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Abstract 



> 

o 

^\i ' The tomographic invertable map of the Wigner function onto the pos- 

^P . itive probabihty distribution function is studied. Alternatives to the 

Schrodinger evolution equation and to the energy level equation writ- 
ten for the positive probability distribution are discussed. Instead of 

f^ ■ the transition probability amplitude (Feynman path integral) a tran- 

sition probability is introduced. A new formulation of the conven- 

C I tional quantum mechanics (without wave function and density ma- 

trix) based on the "probability representation" of quantum states is 

Q^i given. An equation for the propagator in the new formulation of quan- 

tum mechanics is derived. Some paradoxes of quantum mechanics are 
reconsidered. 



Introduction 

During more than 70 years of existence of quantum mechanics, there 
was a dream to reduce misterious and intuitively very unusual notions of 
this theory to the well-known and intuitively acceptable classical notions. 
There was a common prejudice that it is impossible to describe the notion 
of quantum state in the framework of the conventional quantum mechanics 
in terms of the probability density and one is obliged to use either complex 
wave function or density matrix in different representations. Fortunately, 
it turns out that it is possible to associate with a quantum state the usual 



probability density and the use of the wave function or density matrix is not 
mandatory in the conventional quanrum mechanics. This course of lectures 
is devoted to the new formulation of quantum mechanics. 

Quantum mechanics is based on the concept of a complex wave func- 
tion which satisfies the Schrodinger equation |l|]. Several attempts to give 
classical-like interpretations of the wave function were done in p|, |^, Q (see 
also ). These attempts 0, |, Q and related constructions of quasidistribu- 
tion functions in the phase space of the system 0, ^ |], §] give the idea that 
for quantum mechanics it is impossible to describe the state of the quantum 
system in terms of measurable positive probability analogously to the case 
of classical statistical mechanics, where the state of the system is described 
by the positive probability distribution due to the presence of classical fluc- 
tuations. 

Nevertheless, it was shown recently |T^, |Tl], [T2|, |T3|, |1^, |T5[ that in the 



framework of the symplectic tomography scheme [|T^, jT^, which generalizes 
the optical tomography scheme |1^, |T^, it is possible to introduce a classical- 
like description of a quantum state using the measurable positive probability 
(instead of the complex probability amplitude). 

This result was obtained because, in addition to considering a measured 
physical observable in a fixed reference frame in the phase space of the quan- 
tum system, different reference frames in the phase space were considered. 
In the spirit of methodology, it is close to special relativity theory, where 
to get unusual effects due to motion with high velocities, different reference 
frames connected by Lorentz transform must be used. In the quantum case, 
the extra parameters distinguishing different reference frames replace the in- 
formation coded by a phase of the wave function. This approach can be 
considered as introducing a new representation in quantum mechanics which 
can be called the "probability representation" |^0|, ^|. 



The description of the quantum state in terms of positive probability was 
obtained not only for continuous observables like the position [|14|, llB, ^ 



but also for pure quantum observables like spin [^ ^ (see also |^ ). 

A classical formulation of quantum evolution was suggested, and for the 
marginal distribution a new quantum evolution equation was found |T^ which 



is an alternative to the time-dependent Schrodinger equation. This equation 
gives the classical-like description of quantum evolution in terms of a normal- 
ized positive distribution containing complete information on the state of the 
system. Examples of free motion and some excited states of the harmonic os- 



dilator were also considered [|T0], |ri| . The evolution of even and odd coherent 
states |]2S| of a particle in a Paul trap was investigated |2^ in the framework 
of the classical-like description [llO|, . The even and odd coherent states of 



a trapped ion were discussed in recent papers 
states were realized in 



28|. Experimentally these 



29| , |30| . A review of the metod of integrals of motion 



and its application to oscillator's models used in the Paul trap problem is 
given in Ref. pl| . 

the 



The aim of the course of lectures is to discuss, following |13, 14, 15, 32 



notion of a quantum state in the new formulation of quantum mechanics. We 
review the classical-like description of transition probabilities between sta- 
tionary states (energy levels) of quantum systems and obtain analogs of the 
orthogonality and the completeness relations. We show that, if the evolution 
equation describing the dynamics of a quantum system is determined by the 
imaginary part of the system's potential energy considered as a function of 
a complex coordinate, the energy states of the system are determined by the 
real part of the potential energy. The energy levels of the harmonic oscillator 
are rederived in the framework of classical-like alternatives to Schrodinger 
evolution and stationary equations. A new type of eigenvalue problem is 
formulated for the positive and normalized marginal distributions. 



Classical Statistical Mechanics and Tomogra- 
phy Map 

In quantum mechanics, the tomography methods of measuring quantum 
states |T^, |T8|, |19| gave the possibility to introduce new approach to the 
notion of quantum states |^0|, 0, |T^. It turned out that the tomography 
methods can be used in classical statistical mechanics [^, ^]. Following 
32| we start from introducing the positive probability distribution function 
for a state of a classical system. In the course of lectures, we will show that 
both classical and quantum states are described by the same probability 
distribution function (called the marginal distribution function). 

Main expressions for marginal distributions in the optical tomography 
method [0, |19|] as well as in the symplectic tomography method |jl6| are based 



on a theorem which connects the characteristic function with the probability 
distribution function. This connection is valid for quantum states described 



by a density matrix |Q. It is obvious that the same connection also exists for 



classical systems in the framework of classical statistical mechanics. One can 
prove that the Fourier transform of the characteristic function (calculated 
by means of a classical probability distribution) is a positive distribution 
function. We illustrate this statement by an example of a one-dimensional 
system. 

"States" in classical statistics are described by the function / (g, p), which 
is the probability distribution function in the phase space, i.e., 

f{q,p)>0, J f {q, P) dp = P{q) , / / (^' P) ^^ = ^(p) ' 

with P{q) and P{p) probability distributions for position and momentum, 
respectively, {marginals). 

Let the nonnegative function / (g, p) be a distribution function of the 
classical system in the phase space. The coordinates — oo < g < oo and 
— oo < p < cxD are the position and momentum of the system, respectively. 
The function / (g, p) is taken to be normalized 

j f{q,p) dqdp=l. (1) 

We consider an observable X (g, p) which is a function on the phase space of 
the system under study. For the case of classical statistical mechanics, the 
characteristic function for the observable X (g, p) 

X (k) = (e^'^^) (2) 

is given by the relation 

x{k) = Je''''^'^'P^f{q,p)dqdp. (3) 

The Fourier transform of the characteristic function 

wiX) = ^ fx{k)e-^^''dk (4) 
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is a real nonnegative function which is normalized 

fw{X)dX = l. (5) 



In fact, due to Fourier representation of Dirac delta-function, one has 

w{X)=Jfiq,p)S{X{q,p)-X)dqdp. (6) 

The distribution function is nonnegative and the delta-function is also non- 
negative. So we integrate the product of two nonnegative functions over the 
phase space. The result of the integration w {X) is obviously a nonnegative 
function. 

Let us now check the normalization of the function w {X) . We have 

fw{X) dX = f f{q,p) 6{X{q,p)-X) dqdpdX. (7) 

In view of the definition of delta-function, one has 

J5{X{q,p)-X)dX = l. (8) 

This means that 

jw{X)dX = j f{q,p)dqdp. (9) 

Since the distribution function / (g, p) satisfies the normalization condi- 
tion |l]), we have shown that the Fourier transform of the characteristic func- 
tion w {X) given by ^) is normalized too, i.e., it satisfies the normalization 
condition ^. As a classical analog of the quantum symplectic-tomography 



observable introduced in Ref. 16 we consider the classical observable which 



is a linear function on the phase space of the system, 

X{q,p) = M + iyp, (10) 

where the real parameters /i and u are interpreted as the parameters of 
symplectic transform of the position and momentum of the system under 
study. (We discuss only one variable — the position X (g, p) — and do not 
take into account the conjugate momentum.) 

The variable X {q, p) can be considered from two equivalent points of 
view. It can be interpreted as a canonically transformed position which is a 
linear combination of position and momentum in a fixed reference frame in 
the phase space of the system. Another equivalent interpretation of the vari- 
able X {q, p) given by Eq. |10]) is that it is a position of the system measured 



in the rotated and scaled reference frame in the classical phase space of the 
system. 

We use the second interpretation, according to which the real parameters 
/i and V determine the reference frame in the phase space of the system in 
which the position is measured. For the position in the transformed reference 
frame |T0|) , we get from Eq. ^ the distribution function (the tomography map) 



w (X, /x, z.) = -L f ^-^HX-,,-up) f (^^ ^) dqdpdk. (11) 
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Another form for the probability distribution is given by Eq. ^ 

w (X, ^, u) = J f (g, p) 6 {fiq + up-X) dq dp . (12) 

One can see that the marginal distribution is a homogenious function, i.e., 

w (AX, A/i, Xu) = \\\~^w (X, /i, u) . (13) 

We introduced the notation w (X, /i, u) for the probability distribution of 
the position of the classical system in the transformed reference frame in the 
phase space to point out the dependence of the distribution on the param- 
eters /i and u determining the reference frame. Due to the dependence of 
the distribution w (X, /i, u) on these parameters, we call the distribution a 
marginal distribution function. 



The partial case of the canonical transform |TOD is a rotation in the phase 
space 

X = q cos ip + p sin (p . (14) 

This means that we choose the parameters of the symplectic transform 

/x = cos (^ , i/ = sin (yj . (15) 

By introducing the notation for the marginal distribution of the rotated po- 
sition 

w (X, ip) = w (X, fi = cos ip, u = sin ip) , (16) 

we get, in view of [121), 



w (X, (p) = f (g, p) 6 {q cos cp + psin cp — X) dq dp . (17) 



Using Eq. |TID we have another representation for the marginal distribution 
of the rotated position, namely, 

w(X, <^) = — /'e^*^(^~^"°^^-P^''^^)/(g, p) dqdpdk. (18) 
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Introducing the transformed position and momentum 

Q = q cos (p + psin {p; P = —qsin (p + p cos ip , (19) 

in view of the invariance of the volume in the phase space 

dq dp = dQdP, (20) 

we get, using the Fourier representation of delta-function, 

w {X, ip) = 6 {X — Q) f {Q cos ip ~ Psin (f, Q sin cp + P cos (f) dQ dP , 

(21) 



or 



w {X, (p) = f (X cos cp — P sin cp, X sin cp + P cos cp) dP . (22) 



Formula |22|) is mathematically identical to Eq. 12) of Ref. |T9[ where the 
marginal distribution for the homodyne observable was considered. But 
in ^^, the positive classical distribution / (g, p) in the phase space is used 



instead of the Wigner function W (g, p) elaborated in Eq. 12) of Ref. [19 
It is worth noting that the form of the expression for the marginal distri- 
bution w (X, ip) is invariant. The only difference between the quantum and 
classical statistics in the context of the expression for the marginal distri- 
bution w (X, ip) is in the difference between the classical distribution in the 
phase space and the Wigner function. The Wigner function W (g, p) can 
take negative values. The classical distribution function / (g, p) takes only 
nonnegative values. Nevertheless, the result of integration in both cases gives 
the nonnegative marginal distribution w{X, ip) . 
Formula 0) has the inverse 

/ {Qi p) = — / "^ i-^^ /^' ^) ^^P [~^ i^^1 + ^P ~ X)] dX dfi dv . (23) 



In classical statistical mechanics, the admissible marginal distributions in 
formula ^3|) always satisfy the condition that the result of convolution / (g, p) 
is a nonnegative function. 

We have shown that instead of the distribution function / (g, p) the state 
of the classical system in the framework of classical statistical mechanics can 
be determined by the marginal distribution function w {X, fi, z/), in complete 
analogy with the quantum case where the symplectic tomography procedure 



is used 0]. Since the map 



/ {q, p)^w (X, /i, v) 

is invertable, the information contained in the distribution function / (g, p) is 
equivalent to the information contained in the marginal distribution w (X, /i, v) 
For IX = cos (fi and z/ = sin if, we have an analog of the optical tomography 
procedure developed for the quantum case in [|19[. We have to invert for- 
mula 



J2D. The inverse is given by the Radon transform (see, for example, 
Eq. 13) in [^ and also |^ ). For example, if one introduces the distribution 
function in the form 

/ {q, p) = H(i- xo) S{p- po) , (24) 

the marginal distribution takes the form 

w (X, /i, z/) = (5 (X — fixo — upo) (25) 

and 

w (X, (f) = 6 {X — xq cos f — Po sin ip) . (26) 

For classical statistical mechanics, the tomography maps discussed con- 
nect the positive distributions, and in this context our understanding of the 
notion of the classical state for systems with fluctuations is unchanged. 

The evolution equation for the classical distribution function for a particle 
with mass m = 1 and potential U{q), 



df (g, P, t) df (g, p, t) dU{q) df (g, p, t) 
■ +p- 







dt dq dq dp 

can be rewritten in terms of the marginal distribution w (X, /x, z/, t) 



(27) 



d dU ,^ 



d 

"ax" 



0. 



where the argument of the function dU /dq is replaced by the operator 

For the harmonic oscillator with frequency a; = 1, the potential energy term 
U{q) = g^/2 gives in Eq. |28|) the following evolution equation 

dw dw , . 

We used the equality 

1 (90^ f d \'^ d 



and the property 

For the mean value of position in classical statistics, we have 

iX j:/ 



(32) 



{q) = f {q^ p) qdqdp = i w (X, jj, u) e* 6' (jj,) 6 (u) dX djj, dv . (33) 



One can see that in classical statistical mechanics there exists a function 
associated to the position, and by means of this function one can calculate 
the mean value of position using the marginal distribution w (X, /x, v) . 
In classical statistical mechanics, one can introduce the propagator 
H^Y (^5 /^5 ^^ ^' ■• A*'' ^' ■• ^2? ^i) that connects the two marginal distributions 
given for times ti and ti (^2 > ^i) 

w (X, /i, z/, ^2) = / n^^j (X, /i, z/, X', /i', z/', t2, ti) w (X', /i', v\ tx) dX' dfi' du'. 

(34) 
The propagator satisfies the following equation 

= 5{t2-h)5{X-X')5{^-^i')5{v-v'), (35) 



which follows from the evolution equation ^^ . 

Any integral of motion I {q, p, t) in classical statistical mechanics satisfies 
the equation 

dl _ dl dl dH dl dH _ 

dt dt dq dp dp dq ' 

where H (q, p, t) is the Hamiltonian of the classical system. 
Equation ^6|) coincides for 

H = ^^ + U{q) 

with the equation for the classical distribution function |27|) , 

dl (g, P, i) dl {q, P, t) _ dUjq) dl (g, p, t) ^ ^ 
dt dq dq dp 

This follows from the fact that the distribution function itself is the integral 
of motion. 

If one introduces the map 0) for the integrals of motion 

I (X, fi,u,t) = —f e-*^(^-^'?-^P) / (g, p, t) dq dp dk , (38) 

the integral of motion X (X, fi, u, t) satisfies Eq. |2^) in which one has to 
make the replacement w ^ X. 

In classical statistical mechanics, the distribution function / {q, p, t) is a 
function of the integrals of motion, and the propagator that determines the 
evolution of the distribution function has the form 

P {q, P, q', P', t) = S{q'- go (g, P, t)) 6 {p' - po (g, p, t)) , (39) 

where go (g, p, t) and po (g, p, t) are integrals of motion which have the fol- 
lowing property: 

?o(g, P, 0)=g, po{q,p,0)=p. (40) 

Using Eq. ^9|) one can find the propagator for the marginal distribution 
function. 

For example, the initial distribution [2^ ) takes the form 

/o {q, p,t) = 6{p- Po) 6{q-tp- xo) (41) 
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and the initial marginal distribution |25|) reads 

Wo (X, /i, u, t) = 6 {X - fitpo - fixo - upo) 



(42) 



In quantum and classical statistical mechanics, the forms of the propa- 
gators determining the evolution of the marginal distributions w {X, /z, u, t) 
are identical for linear systems like an oscillator or free motion. 



New Notion of Quantum State 

We consider now a new approach to the notion of quantum state. It was 
shown [|l^ that for the generic linear combination of quadratures which is a 
measurable observable {% = 1) 



X = fiq + up . 



(43) 



where q and p are the position and momentum, respectively; the marginal 
distribution w{X, fi, v) (normalized with respect to the variable X), de- 
pending on the two extra real parameters /x and v, is related to the state of 
the quantum system expressed in terms of its Wigner function W{q, p) as 
follows: 

dk dq dp 



w (X, n, u) = / exp [—ik{X — ^q — vp)] W{q, p) 



i2ny 



(44) 



We use the same notation as in the classical case. If one has a pure state 
with the wave function \1/ (y), the marginal distribution has the form found 
in Ref. 



w {X, /i, u) 



1 



27r z/ 



^ {y) exp 



2u 



lyX' 



V 



dy 



(45) 



The physical meaning of the parameters /i and u is that they describe an 
ensemble of rotated and scaled reference frames in which the position X is 
measured. For fi = cos (f and u = sin ip, the marginal distribution ^) is 
the distribution for the homodyne-output variable used in optical tomogra- 
phy ||19[|. Formula ^) can be inverted and the Wigner function of the state 



can be expressed in terms of the marginal distribution |T^ : 

W^(Q') p) = — 1^ (-^! /^) ^) exp [—i {fiq + vp — X)] d^ dv dX . 
27r J 



(46) 
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Since the Wigner function determines completely the quantum state of a 
system and, on the other hand, this function itself is completely determined 
by the marginal distribution, one can understand the notion of the quantum 
state in terms of the classical marginal distribution for squeezed and rotated 
quadrature. 

So, we say that the quantum state is given if the position probability dis- 
tribution w (X, /i, v) in an ensemble of rotated and squeezed reference frames 
in the classical phase space is given. 

It is worth noting, that the information contained in the marginal distri- 
bution w (X, yU, v) is overcomplete. To determine the quantum state com- 
pletely, it is sufficient to give the function for arguments with the constraints 
(/i^ + 1/^ = 1) which corresponds to the optical tomography scheme |1^, |37 



i.e., /i = cos if and the rotation angle 99 labels the reference frame in the 
classical phase space. 

So, we formulate also the notion of quantum states as follows: 
We say that the quantum state is given if the position probability distri- 
bution w {X, if) in an ensemble of rotated reference frames in the classical 
phase space is given. 

Since the quantum state is defined by the position distribution, one could 
associate an entropy with the state using the standard relation known in 
classical probability theory, i.e., the entropy S {fi, v) is given by the formula 

S{^,u) = -fdXw {X, /x, z/) In [w (X, //, z/)] . (47) 

If we use the distribution w {X, ip) , the entropy S {ip) depends only on the 
rotation angle. 

The discription of quantum states by the probability function gives the 
possibility to formulate quantum mechanics without using the wave function 
or density matrix. These ingredients of the quantum theory can be considered 
as objects which are not mandatory ones since they are not directly measur- 
able. The marginal probability distribution function w (X, /i, z/), which can 
be measured directly, replaces the wave function in the new formulation of 
quantum mechanics. Since the quantum mechanics formalism is reduced to 
the formalism of classical probability theory, well-known results of the prob- 
ability theory can be used to get new results in quantum theory (including 
quantum computing, teleportation, and quantum cryptography). 
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One can also use the introduced formulation of the notion of quantum 
states to describe situations in which the states are either close or essen- 
tially different. We say that two states are close if their distributions are 
close, i.e., all the highest momenta of the distributions differ very slightly. 
We also say that two states are substantially different if their distributions 
differ substantially, i.e., there are highest momenta for the two distributions 
with large corresponding differences. The notion of distance in quantum 



mechanics using the tomography map was discussed in [38 



Quantum Evolution and Energy Levels 



As was shown in [0, for systems with the Hamiltonian 

S2 



H 



P 



+ V{q) 



(48) 



the marginal distribution satisfies the quantum time-evolution equation, be- 
ing the integral equation determined by the imaginary part of the potential 
energy considered as a function of a complex coordinate. The evolution 
equation reads 



d 

w — ib—-W 
Oh' 



V 



V 
1 



1 d 



V d 

I — - 



d/dXdfi 2dX) 
w = 0. 



d V d 



d/dXd^i 2dX 
This equation is alternative to the Schrodinger equation 

i^ = H^H 
and it can be obtained from the equation for density matrix 

p + i[H, p]= 0, 
in view of the following formulas: 



qW{q,p) 




d_ 
dp 



w (X, /i, u) 



(49) 



(50) 



(51) 
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and 



d d 

Q-W{q,p) > // — w (X, /i, I/) , 



f d Y' d 
pW{q,p) — > ~[dx) ^^(^'/^'^) 



-p(X,X') -^ U7^-^H^(^'^) 



(52) 



aX'^' ' ' \2 dq 



xp(x, X') -^ fg + il-1 w^(g,p) 



(53) 



2 9j9^ 

x'p(x,x') ^ (^-^I)^(^'P)- 

Equation 49) can be considered as a Fokker-Planck-like equation of classi- 



cal probability theory. The measurable position is a cyclic variable for the 
evolution equation. 

In order to compare the classical and quantum evolution equations, let 



us rewrite the quantum evolution equation ^Q]) in the form of a series, 



Here 

f)2n+lY 

^""' (^) = ^^^ (^ = ^) ' (55) 

where the operator q is given by Eq. p9|). 

Equation ^5]) can also be presented in the form 

„_,^_^ ,^) 9 „+2 V r^ f ^4V"" (-ir ■- = 0. (56) 
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The three first terms give the ^ — ;> classical Boltzman equation. It is 
important that both classical and quantum evolution equations are written 
for the same function w (X, /i, z/). 

Let us rewrite, following |jl3|, the Schrodinger equation for the stationary- 



state density matrix pE of the quantum system with Hamiltonian ^81) 

HpE = PeH = EpE (57) 

in terms of the time-independent marginal distribution we {X, p, v) of the 



squeezed and rotated quadrature introduced in [16|. We have 



+Re1/ 



AdX, 
d 



d^ 



r^x 



9z/2 
dX) dp 



We- -p 



d^ 



We 



= Ewe ■ 



(5^ 



The positive marginal distribution (eigendistribution) satisfies this eigenvalue 
equation and also the equation 



d 

fJ'TrWE 
Oh' 



2lmV 



i d 
2^dX 



dX 



d_ 
dp 



We- 



(59) 



Equation |59|) follows from the evolution equation |49|) for the marginal dis- 
tribution of the quantum system (see Ref. |TD[ ), if the marginal distribution 
does not depend on time. Thus, the normalized marginal distributions of 
stationary states of quantum systems satisfy the system of two equations |S3) 



and p9|) . 

We consider an example of the quantum harmonic oscillator since it is one 
of the most important quantum systems. For this case, using the Hamiltonian 



2 2 



we reduce Eq. |^) to the following one (in view of Ref. [10] ): 



d d 

w — p-—w -\- V -— w = 
av dp 



(60) 



(61) 
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The marginal distribution of the oscillator's ground state is 



,(os) 



Wq {X, /i, u) 



1 



TT (/i2 + U^) 



exp 



X' 



fi' + u' 



(62) 



The marginal distribution must be consistent with the uncertainty relation, 
i.e., 



fw{X, 1, 0)X'^dX~nw{X, 1, 0)XdX 
fw{X, 0, l)X'^dX-Uw{X, 0, l)Xdx\ 



>4. (63) 



Propagator 



In Ref. Ill 



the classical transition-probability density from an initial 
position X' measured at time t = in the reference frame in the classical 
phase space labeled by the parameters //'; z/' to the position X measured 
at time t in the reference frame in the classical phase space labeled by the 
parameters /i; z/ was introduced. This classical transition-probability density 



is the propagator for the evolution equation ^9D for the marginal distribution 
and the propagator is the kernel of the integral relation 

w (X, fi,u,t)=fu (X, /i, z/, X', fi', u', t) w {X', fi', z/', 0) dX' dfi' du' . 

(64) 
The classical propagator has a specific feature, it takes into account that 
the transition probability is considered in different references frames in the 
phase space. In view of this fact, parameters of reference frames fi and z/ are 
present in the evolution equation. Due to this, the equation for the propaga- 
tor slightly differs from the Smoluchowski-Chapman-Kolmogorov equation 
elaborated in the classical probability theory. 

The classical propagator can be related to a quantum propagator (the 
Green function) for the density matrix p {X, X', t) in the coordinate repre- 
sentation. For a pure state with the wave function \1/ (X, t) , we have 



p(X,X', t)=v[/(X,t)vl/*(X', t). 



(65) 
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The Green function of the Schrodinger equation G (X, X', t) connects the 
wave functions at the initial time moment t = and at time t 



m{x,t)= I G (X, X', t) m (X') dx' . 



(66) 



We have for the density matrix ^), in view of relation ^6|), the following 
expression: 



p (X, X', t)= j K (X, X', r, Y\ t) p {Y, y, t = 0)dY dY' 



(67) 



where the propagator K (X, X', Y, F' , t) is expressed in terms of the Green 
function (for unitary evolution) 



K (X, X', Y, Y', t)=G (X, r, t) G*(X', Y' ,t) . 



(6^ 



Since the relation of the density matrix to the marginal distribution is known 
for any time t (given by |65|) and |66|) ), it is possible to obtain 



K{X, X\ Z, Z\ t) = J^,jy^ exp |. (r-/i 

— % y + « u > 

v' 2z/' ^J 

X n (F, /i, X - X', y, ^', z/', t) d/i d^' dY dY' dv' . 



X + X'^ 



(69) 



Thus, given the classical propagator for the classical marginal distribution, 
the propagator for the density matrix is also given. Formula ^9|) can be 
converted. 



Deriving formula |69|) we used the relations 
^ (<?; p) = t: — / "^ (-^? /^5 ^) exp [—i {pq + up — X)] dp du dX . 

2 TT J 



p(X,X') = ^ fwiY,p,X-X')exp 
Zn J 



I \Y~p 



X + X'^ 



dpdY , 



and 



w (X, /i, v) 



r z/ J 



27r, 
X exp 



—2 X — U 



dZdZ'. (70) 
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The last formulas give some relationships between the marginal distribution 
w {X, /i, z/), the Wigner function, and the density matrix in the coordinate 
representation. 

In Ref. @, the Wigner function was introduced in terms of the density 
matrix 

W{q, p) = /P (9 + ^' 9- I) ^"'"'du, (71) 

which can be rewritten as 

Wiq,p) = f p{Z, Z') 6(z-q--) 6(z' -q + -]e~'P''dudZdZ', 



{q,p)=2 I p (Z, Z') e-2^P(^-'^) 5 {Z' + Z - 2q) dZ dZ' . (72) 



or 

W 

Comparing formulas |72D and |70|) one can conclude that the Wigner qua- 
sidistribution function W (g, p) and the classical probabihty distribution 
w (X, yU, v) , the latter being a positive and normalized function, are ob- 
tained using similar integral transforms of the density matrix. 

The difference between the two functions is determined by the difference 
in the kernels of the integral transforms. In the case of the Wigner transform, 
the kernel reads 

Kw {Z, Z', q,p)=2 e-2^P(^-'?) 6 {Z' + Z - 2g) . (73) 



In the case of the symplectic tomography transform suggested in Ref. [16 
the kernel reads 



Km {Z, Z', X, p, v) = exp 



.Z-Z' / Z + Z'' 



■ (74) 



Due to the difference of the kernels, the Wigner function takes negative values 
and the marginal probability distribution is nonnegative function. 

If one writes the classical propagator as a function of the initial time 
moment ti and the final time moment ^2 (i-e., ti 7^ 0), relation ^4]) can be 
rewritten as 

w (X, /i, z/, ^2) = / n (X, yU, u, X', /x', z/', ^2, h) w (X', /x', z/', ti) dX' dp! du'. 

(75) 
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From the physical meaning of the classical propagator, the nonlinear integral 
relation follows 

n (X, /i, u, X', /i', u', t2, ^i) = / n (X, /i, u, X", /x", u", h, t') 

X n (X", /i", z/", X', /i', z/', t', ti) dX" dfi" du" . (76) 

This relation means that if the system is initially located at the point X' at 
time ti in the reference frame in the phase space labeled by the parameters 
/i'; z/', the probability for the system to arrive at the point X in the reference 
frame in the phase space labeled by the parameters /z; u at time ^2 is equal to 
the probabilities to arrive at an intermediate point X" in the reference frame 
in the phase space labeled by the parameters /i"; z/" at time t' integrated over 
all the intermediate positions and all the intermediate reference frames. 

The above integral equation |76D is an analog of the Smoluchowski-Chapman- 
Kolmogorov relation which in the approach introduced in Ref. [|l^ is gener- 



alized to the case of families of conditional probabilities if different reference 
frames in the phase space (parameters /i and u) are taken into account. Also 
the propagator satisfies the differential equation (see Ref. [jlT| ) 

dt2 diy 



\ d/dXdfl ^2dX^ 



( 1 d_ .u_d_' 

[ d/dXdJ1^^2dX 



n 

= 6{t2- h) 6{X- X') 6{fi- fi') 6{u- v') . (77) 

The classical propagator satisfies the initial condition 

n (X, fi, z/, X', /i', z/', t,t) = 6{X- X') 5 (/i - /i') 6{u- v') . (78) 

The relation that could be used to express the classical propagator in 
terms of the functional integral can be also written 



n (X, /i, z/, X', /i', z/', tf , tj^) 

~ / 11 {n(Xfc_|_i, /ifc+i, z^fc+i, X]^, /Xfc, z/fc, tfc_|_i, tk) dX^dfikduk} (79) 
•' k=i 

where the time interval t^ — tj^ = N t ; t^ = t^^ + kr ; k = 1, 2, . . . , N. 
Taking in relation |7^) the limit r -^ ; X — > cx), one obtains the expression 
for the classical propagator in terms of the functional integral. 
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Quantum Transition Probabilities 

We express quantum-transition probabilities in terms of an overlap in- 
tegral of classical-like marginal distributions describing the initial and final 
quantum states. 

If the initial pure state of a quantum system is described by the marginal 
distribution 

Win = ^1 (^' /"' '^) (80) 

and the final state of the quantum system is described by the marginal dis- 
tribution 

Wf = W2 (X, fi, v) , (81) 

the probability of the quantum transition P12 (1 =^ 2) can be obtained using 
the known expression for the probability in terms of an overlap integral of 
the Wigner functions Wi [q, p) and W2 {q, p) of the initial and final states 
(see, for example, Ref. [^ ) 

Pi2 = 7^ fWiiq,p)W2{q,p)dqdp. (82) 

Zn J 

For the transition probability, one can obtain, in view of relation ^), the 
expression in terms of the marginal distributions 

P12 = / wi{X, /i, z/) W2{Y, -/i, -z/) exp [i (X + Y)] ^^— . (83) 



As follows from relation ^), any pure normalized quantum state is described 
by the marginal distribution wp {X, fi, v) , which satisfies the additional con- 
dition 

/■ /^^ ^ ,^ r NT /,^ ^^.^ dududX dY , ,, 

y «;p(X, /i, v) wp{Y, -/i, -v) exp \i (X + F)] -^^^- = 1 . (84) 

The complex wave functions, which belong to different energy levels of a 
quantum system, are orthogonal. This orthogonality condition is expressed 
in terms of the classical marginal distribution as the relation 

/■ /^A \ /T^ \ r-/^A ^^.^d^idvdXdY ^ , ^, 

/ Wn\X, /i, v) Wm\Y-, -fi, -v) exp [i (X + Y)\ = dmn , (85) 
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where the labels m, n correspond to the energy levels Em, En. The pure 
states I n) satisfy the completeness relation 

Y, \n){n\=l. (86) 

n 

This relation can be rewritten as the condition for the marginal distributions 
of the pure states with the energies En 

5:«;„(X, /i, z/)=^(™)(X,/i, z/), (87) 

n 

where the distribution w^^^^ describes the white noise, 



(wn) r V .. ..\ I dx dy dk 



w^^''> (X, /i, z/) = ^^^V^ "^P 

Z TV 



ik (/ix + uy) — ikx — ik"^ — 



Thus, the classical marginal distributions describing the energy levels 



of quantum systems are positive solutions to the system of equations |5^) 
and |SD|) and these solutions satisfy the orthogonality condition ^) and the 
analog of the completeness relation ^). The distributions form an interesting 
mathematical set that differs substantially from the usual Hilbert space of 
states described by the normalized complex wave functions. Of course, the 
structure of the Hilbert space can be traced using the map, which connects 
the states expressed in terms of the density matrix and the states expressed 
in terms of the marginal distribution functions. 



Propagator for Systems with Quadratic Hamil- 
tonians 

As an example, we consider the system with the quadratic Hermitian 
Hamiltonian 

H = ^iQBQ)+CQ, (89) 

where one has the vector-operator Q = {p, q) . The symmetric 2x2 matrix 
B and real 2-vector C depend on time. The system has linear integrals of 



motion (see Ref. 11, 11 ): 

/ (t) = A (t) Q + A (t) . (90) 
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Here the real symplectic 2x2 matrix A (t) and the real vector A (t) satisfy 
the equations 

A = iABay, A = iAayC, (91) 

and the initial conditions 

A(0) = 1; A(0) = 0. (92) 



As follows from relation |J) and from the property of the Wigner function 



(see Ref . |^ ) , the classical propagator is 

n (X, fi, u, X', fi', u', t)=6{X-X' + AT A^^A) 6 {M' - U k'^) , (93) 
where the vectors A/" and A/"' are 

^f = {v, /i) , 

M' = (^',/i'). 

For the quadratic systems without linear terms (C = 0) , the classical prop- 
agator is 

n (X, /i, u, X', fi', u', t) = S{X- X') S {AT' -ATA-^) . (94) 

Thus, if one knows the linear integrals of motion, i.e., the matrix A (t) and 
the vector A (t), one knows the classical propagator. 
For free motion with the Hamiltonian 

H = ^, (95) 



there are two linear invariants found in Ref. [|39| , ^ 

Po{t)=P, qoit)=q-pt. (96) 

This means that A (t) = 0, and the symplectic 2x2 matrix reads 

Thus, we have 

Af A-\t) = {u + fit, fi) . (98) 

22 



Consequently, the classical propagator of free motion has the form 

U^^\X, /i, z/, X', /i', z/', t) = 6{X- X') 6{u' -u- /it) 6 {12- /i') . (99) 



For the harmonic oscillator with the Hamiltonian ^0]), linear invariants 
are known (see Ref. ^^ ^ ), and the matrix A (t) is 



A(t) 



cos t sin t 
— sin t cos t 



(100) 



This means that for the harmonic oscillator 

J\f A^^{t) = (z/ cos t — /i sin t, z/ sin t + /i cos t) . (101) 

Consequently, the classical propagator of the harmonic oscillator is 

n(os)(x,/i, z/,x',/i', z/') 

= 6 {X — X') 5 {v — V cos t + /i sin t) S {fi' — u sin t — fi cos t) (102) 



Energy Levels of the Harmonic Oscillator 

The marginal distribution of the coherent state of the harmonic oscillator 
has the form obtained in Ref. |TT| : 



WaiX, H, v) 



[vr (/x^ + u')] 



2N1-1/2 



exp 



— \a\ 



X2 a''{v + i^f 



/i2 + zy2 2(/i2 + Z/2) 



a*2(z/-i/i)2 i\/2aX{i' + ill) iy/2a*X{u - ifi\ 



2 (/i2 + zy2 



/i2 + Z/2 



/i2 + Z/2 



03) 



The eigendistribution function for the energy level of the harmonic oscillator 
satisfies the eigenvalue equation 



dv 



'd_' 



MT-l^^'-^'^iM' 



Ewe{X, fi, z/) 



WEiX, jj,, v) 

(104) 



This equation can be rewritten for the Fourier component of the marginal 
distribution 



WE(k, /i, z/) = — / we{X, ij,, v) exp {—ikX) dX 
zn J 



(105) 
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in the form 

1 



2F 



dv 



+ 






k {^ + I' ) > Wsik, fi,iy) = E wsik^ /i, u) 



(106) 

Since the marginal distribution of the stationary state of the harmonic oscil- 
lator must satisfy the stationarity condition found in Ref. |T^ , 



d d\ 

f^^r.~^g-j we{X, /x, z/) = , 



du 



(107) 



Eq. p,06| ) is equivalent to the equation for axially symmetric wave functions 
of a two-mode harmonic oscillator with mass m = k"^, frequency u = 1/2, 
and angular momentum M = 0. The wave function corresponding to zero 
angular momentum is expressed in terms of the Laguerre polynomials 



Wn{k, /i, z/) = -— exp 

Z 71 



I^\li' + P' 



2, ,2' 



Lr 



/C yU + k u 



(108) 



The main quantum number n of the one- dimensional harmonic oscillator 
under discussion is equal to the integer radial quantum number n^ of the 
artificial two-mode oscillator 



n = Ur 



rir 



0,1,2,.... 



(109) 



The energy level of the artificial symmetric two-mode oscillator labeled by 
the radial quantum number n^ and the angular momentum M as 



E, 



nr,M 



u{\M\ + l + 2n^) 



(110) 



for uj = 1/2 , M = , rir = n gives exactly the spectrum of the one- 
dimensional oscillator 



E„ = E. 



rir , M 



n + 



1 



n = 0, 1, 2, 



(111) 



To find the marginal distribution, we have to calculate 



Wn{X, fi, Z/) 



1 

2^ 



exp 



P(^2 ^ ^2 



+ ikX 



L.(^^^),k. 
(112) 
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In view of the integral 

— / exp h ikX \ Ln\ — \ dk 

= vr-i/2 2-" {n\)-' exp {-X') HI{X) , (113) 

one obtains the marginal distribution 

u.„(A', fi, (.) = lir{i/ + i,'-)\-"'2-(n\)-' 

It is worth noting that the normalization condition for the marginal distri- 
bution Wn{X, /i, v) implies the condition for the Fourier component 

/ Wn{k-, /i, v) exp (ikX) dkdX = 2nWn{k = 0, fi, u) = 1 . (115) 



We take solutions |108|) without using the normalization condition in terms 
of the variables /i and u of the artificial two-mode oscillator, but using the 
normalization condition of the marginal distribution in terms of the variable 
X and the corresponding property of its Fourier component. 



Tomography of Spin States 

Let us introduce the probability distribution for the spin projection in a 
given direction considered in a rotated reference frame. For arbitrary values 
of spin, let the spin state have the density matrix 

Pmm' = U^ I P^^^ I J^') > m = -j,-j + 1,...J -IJ, (116) 



where 



J3 I jm) = m I jm) , 

f I jm) = jij + 1) I jm) 



(117) 



and 



P^'^= E E P^£m'\jm){jm'\ . (118) 

m=—j m'=—j 
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The operator p'^^'> is the density operator of the state under discussion. The 
diagonal elements of the density matrix determine the positive probability 
distribution 

P^L = woim) , (119) 

which is normalized, 



J2 ^o( 



771} 



1. 



(120) 



m=-j 



In Refs. [^ ^, ^, a general group construction of tomographic schemes was 



discussed, and this scheme was also used for spin tomography in Refs. [22, 23 



The idea is to consider the diagonal elements of the density matrix in another 
reference frame. The density matrix in another reference frame reads 



P\ 



U) 



mi m2 



\ / mi 



:i2ii 



m2 



Here the unitary rotation transform T> depends on the Euler angles a, /?, 7 
and, by definition, the diagonal matrix elements of the density matrix yield 
the positive normalized probability distribution. For the diagonal elements 
of the density matrix |121|) , 

mi = 7712 ■ 



We introduce new notation and rewrite equality |121| ) for rrii = 7712 in the 
form 



w I 



(i) 



i(i)* 



mi,a,/3,7)= E E ^mU("'^'^) /^™U ^-'U^"'^'^) " ^^22) 



m-^=-] 171^=-] 



Here the matrix elements D^ ^, (a, /3, 7) (the Wigner function) are the ma- 



trix elements of the rotation-group representation 



i(i) 



D'^'m{c^.M)=e''-^di:,^{(3)e^ 



jU) 



(123) 



where 



"-m'miy) 



1/2 



(j + 7n')\{j — 7n')\ 
{j + 7n) ! (j — 7n) ! 

cos P) 



cos 



/3^ 



m'+m 



sm 



f3' 



j-)(rn' —m,m' +m) 
j—m' 



(124) 



26 



and P^"'^)(x) is the Jacobi polynomial. 
Since 

the marginal distribution depends only on two angles, a and /3. 
Thus, let us denote 

w (mi, a,/?) = w (mi, a,/?, 7) , (126) 

which satisfies 

^ w(mi,a,/5) = l. (127) 

mi=-j 

For a spin-1/2 state with spin projection +1/2 and wave function 



V'+i/ 



2 



or with density matrix 

1 



^+" VO 
the marginal distribution is equal to 

w(-,a,pj =008"^ - for mi = +- , (128) 

and, correspondingly, 

w ( — -, a,/? ) = sin^ — for mi = — - . (129) 



2' " / 2 2 

In Ref . PH] , by using the properties of the Wigner function and the Clebsch- 



Gordan coefficients, formula |122|) was inverted and the density matrix was 



expressed in terms of the marginal distribution 

2j J3 

E E 

J3=0 m3=-J3 mi=-j' 



;-l)-'. ^ ^ (2J3 + 1)' E y(-l)'"^^(mi,a,/5) Z^E(«,A7) 



X ( ^' ^' -^'^^ V ^' •^' 3^\^ = pU) ^ n 30) 

\m,i — ?7ii / V'^i """^2 ^^3/ Svr^ f^m'^m'^- 
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Here m,m = —j, —j + l,...,j and one integrates over rotation parameters 
a, f3, 7. 

To derive formula |13CI|) , we used the known property of the Wigner func- 
tion: 



D 



O'l) t, ,\ n(h) I. ,\ nUa) / ^ duj _ f j^ j^ jg 



'^U^)D'^U^^^^rnM-' 



r2 



',n^ \m[ m'2 TTig 



^ ( Jl J2 J3 ) ^^3^) 

* nil ^TT'2 ^^3 



where 



f dto =[ "^ da r sin p dp [ "^ dj, (132) 

J Jo Jo Jo 

along with the orthogonality property of 3j-symbols: 



V- ^^ f Ji 32 ] \ f ]l ]2 ]' 



^ -^ ' "^ . '^ . \mi nio -mj \mi nio -m' J " ' 

'E E (2J + 1) ( ^' ^' ^ ) ( ^\ ^' ^ ) = 5m,m', S^.m',- 

^^ , '^^ . \mi mo —ml \m\ rrin —m, J '- ^ ^ ^ 

Formula p.3CI| ), being the inverse of p.22| ), is an analog of the Radon transform 
for spin states. Given a measurable marginal distribution for arbitrary spin, 
one can reconstruct the state density matrix by means of this relation. 

The results obtained enable one to measure the spin state by measuring 
a spin projection on a given axis. One obtains the experimental probability- 
distribution function which depends on two angles determining the axis. Us- 
ing the relationship between the probability distribution and the state density 
matrix, one reconstructs all the information about the quantum spin state. 
This means that the probability distribution can be used instead of complex 
spinors and density matrices for the spin-state description since it contains 
complete (even overcomplete) information on the state. 



Quantum Measurements and Collapse of Wave 
Function 

We will review the discussion of quantum measurements done in Refs. [ffO, 



Tl| , |14| . It is known (see, for example, Refs. |^4|, ^] ) that quantum mechanics 
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is problematic in the sense that it is incomplete and needs the notion of a 
classical device measuring quantum observables as an important ingredient 
of the theory. Due to this, one accepts that there exist two worlds: the 
classical one and the quantum one. In the classical world, the measurements 
of classical observables are produced by classical devices. In the framework 
of standard theory, in the quantum world the measurements of quantum 
observables are produced by classical devices, too. Due to this, the theory of 
quantum measurements is considered as something very specifically different 
from classical measurements. 

It is psycologically accepted that to understand the physical meaning of 
a measurement in the classical world is much easier than to understand the 
physical meaning of analogous measurement in the quantum world. 

As was pointed out in Refs. |Ty, |TT|, all the roots of the difficulties of 



quantum measurements are present in classical measurements, as well. Using 
the relations of the quantum states in the standard representation and in the 
classical one (described by classical distributions), one can conclude that 
complete information on a quantum state is obtained from purely classical 
measurements of the position of a particle made by classical devices in each 
reference frame of an ensemble of classical reference frames, which are scaled 
and rotated in the classical phase space. 

These measurements do not need any quantum language if we know how 
to produce, in the classical world (using the notion of classical position and 
momentum), reference frames in the classical phase space differing from each 
other by rotation and scaling of the axis of the reference frame and how to 
measure only the position of the particle from the viewpoint of these dif- 
ferent reference frames. So, knowing how to obtain the classical marginal 
distribution function w (X, /i, z/) which depends on the parameters /x; v, la- 
beling each reference frame in the classical phase space, we reconstruct the 
quantum density operator. 

Thus, we avoid the paradox of the quantum world which requires for its 
explanation measurements by a classical apparatus accepted in the frame- 
work of standard treatment of quantum mechanics. But the difficulties of 
the quantum approach are present, since we need to understand better the 
procedure of measurement in a rotated reference frame in the phase space of 
the classical system. The problem of wave function collapse Q, ^ reduces 
to the problem of a reduction of the probability distribution which occurs 
as soon as we "pick" a classical value of the classical random observable in 
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the classical framework of [0, |Tl|. This means that we "solved" the para- 
dox of the wave function collapse reducing it to the problem of standard 
measurement of a classical random variable used in the probability theory. 



The approach developed in |T^, |TT| enables one to transform such an un- 



pleasant problem of standard quantum mechanics as the need of a classical 
device and the reduction of wave packets into the standard problem of classi- 
cal measurements of classical random variables. In fact, this means that the 
problem of classical measurements is as difficult as the problem of quantum 
measurements. An important analogy with methodology of special relativ- 
ity arises: It turns out that it is necessary to introduce a consideration of 
events in the set of moving reference frames in space-time in order to ex- 
plain relativistic effects, and it is necessary to introduce a consideration of 
events in the set of rotated and scaled reference frames in the phase space 
in order to explain the nonrelativistic quantum mechanics in terms of only 
classical concepts of classical fluctuation theory. But these reference frames 
are the reference frames in the phase space (not in space-time). Possibly, a 
combination of these two approaches can be generalized to give a classical 
description of relativistic quantum mechanics. 

One can conclude that the stationary states of quantum systems (for 
example, of a harmonic oscillator) can be obtained using classical-like alter- 
native equations to the Schrodinger equations. A new type of eigenvalue 
problems for real positive marginal distributions is formulated. The analogs 
of orthogonality and completeness relations for the wave functions are for- 
mulated in terms of conditions for the marginal distributions as well as the 
transition probabilities among the energy levels. The criterion for determin- 
ing the pure states of the quantum system is given in terms of the classical 
marginal distribution. 

Thus, using the marginal distribution one can formulate the standard 
quantum mechanics without the complex wave function and density matrix. 
But the position distributions in an ensemble of classical reference frames in 
the phase space play an important role. 

It should be pointed out that in the standard formulation of quantum 
mechanics there exist different representations such as the coordinate repre- 
sentation, the momentum representation, etc. The counterpart of this variety 
in the classical formulation is related to different tomography schemes like 
optical tomography |]TP|, |^, symplectic tomography [^ |TB[, and photon 



number tomography [^, ^ ^. The photon number tomography uses the 
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marginal distribution of a discrete variable, which corresponds to a number 
representation. Just as different representations in the standard formulation 
of quantum mechanics are related by some transformations in the Hilbert 
space of states, the marginal distributions of different tomography schemes 
can be transformed into each other. This transformation consists of two 
steps. First, one makes a map of the marginal distribution (in one of the 
tomography schemes) onto the Wigner function and then one makes another 
map (of the different tomography scheme) of this Wigner function onto the 
corresponding marginal distribution. 

The construction introduced in spirit is similar to the Moyal approach [^ 



which considers quantum mechanics as a statistical theory. But in Ref. | i5[] , 
the quantum state was described by a quasidistribution function in the phase 
space that is identical to the Wigner function. Thus, the "negative probabil- 
ities" to find the system in some domain in the phase space is an unavoid- 
able feature of the Moyal approach. The density matrix was introduced in 
Ref. [^, |5^ . In the framework of the new formulation of quantum mechanics, 
the density matrix is not mandatory to be used. 

In the introduced formulation of quantum mechanics, only positive prob- 
abilities of the measurable position in an ensemble of reference frames in 
the phase space of the system is used. It is remarkable that in the positive 
probability representation the states in quantum mechanics are described 
identically with the states in classical statistical mechanics if one uses the 
positive marginal distribution w (X, /i, z/) (though the sets of the distribu- 
tions in the classical and quantum cases are different). 

The difference between classical statistical mechanics and quantum me- 
chanics in the formulation introduced appears in the dynamics of the marginal 
distributions, since in quantum mechanics the evolution equation for the pos- 
itive probability distribution has a different form from that in the classical 
case. It is remarkable that the relations of the propagators (conditional 
probabilities) in the phase space representation and in the probability rep- 
resentation are described by the same formula both in classical statistical 
mechanics and in quantum mechanics. 

For linear systems (oscillators), the propagators in classical stastistical 
mechanics and in quantum mechanics coincide. The difference for these sys- 
tems in the quantum and classical cases is due to the fact that not all positive 
probability distributions w (X, /i, v) are realized for classical systems. Also 
not all marginal distributions are admissible in the quantum case, but only 
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which satisfy uncertainty relations. 

We have demonstrated that spin states and states of a trapped ion ||^ can 
be described by measurable positive probability distributions. This implies 
that quantum-mechanical systems can be considered in the framework of the 
same formalism of probability theory as classical statistical systems. Thus, 
the known results of classical probability theory can be applied to the study 
of quantum states. For example, the central limit theorem can be used for 
describing multimode systems. The approach developed can be elaborated 
for solving many problems of quantum optics [^ |S2|, ^ and quantum com- 
puting. It is important to study the Schrodinger uncertainty relation |^4[ in 
the framework of the new approach. Linear integrals of motion for quadratic 
systems [|55| , |56| are useful to obtain the propagator of the new evolution 
equation for the marginal distribution |E^. The new approach can be also 
applied to study nonlinear coherent states [0, ^, |6^ . 

We have shown that quantum mechanics can be formulated without wave 
function and density matrix using the tomographic probability representa- 
tion. The general approach to the tomographic map and relations among 
different tomography schemes are discussed in Refs. [^ 
phy of spin states for two particles is described in Ref. 
new approach to quantum mechanics is given in Ref. |]B^ 



5^]. The tomogra- 

63[ . A review of the 

The generalization 

One 



65|. 



of the metod to the case of the field theory one can find in Ref. 

can conclude that the problem of formulation of quantum theory using only 

probabilities both for continuous and discrete observables has the solution 

in the framework of the tomographic probability representation of quantum 

mechanics. 



Acknowledgments 

The author thanks the Organizers of the XXXI Latin-American School 
of Physics for invitation to give a course of lectures and El Colegio Nacional 
for kind hospitality. 

This study has been partially supported by the Russian Foundation for 
Basic Research under Project No. 17222. 



References 



32 



[1] Schrodinger, E., Ann. d. Physik (Leipzig) 79, 489 (1926). 

[2] De Broglie, L., Compt. rend. 183, 447 (1926); 184, 273 (1927); 185, 
380 (1927). 

[3] Madelung, E., Zeits. f. Physik 40, 332 (1926). 

[4] Bohm, D., Phys. Rev. 85, 166; 180 (1952). 

[5] Terletskii, Ya. A., and Gusev, A. A., (Eds.), Problems of Causality in 
Quantum Mechanics (Inostrannaya Literatura, Moscow, 1959) [in Rus- 
sian] . 

[6] Wigner, E., Phys. Rev. 40, 749 (1932). 

[7] Husimi, K., Proc. Phys. Math. Sac. Jpn 23, 264 (1940). 

[8] Glauber, R. J., Phys. Rev. Lett. 10, 84 (1963). 

[9] Sudarshan, E. C. G., Phys. Rev. Lett. 10, 277 (1963). 

[10] Mancini, S., Man'ko, V. I., and Tombesi, R, Phys. Lett. A 213, 1 (1996). 

[11] Mancini, S., Man'ko, V.. I., and Tombesi, P., Found. Phys. 27, 801 
(1997). 

[12] Man'ko, V. I., "Energy levels of a harmonic oscillator in the classical- 
like formulation of quantum mechanics," in Dremin, I. M., and 
Semikhatov, A. M., (Eds.), Proceedings of the Second International 
A. D. Sakharov Conference on Physics, Moscow, May 1996 (World Sci- 
entific, Singapore, 1997), p. 486; "Optical symplectic tomography and 
classical probability instead of wave function in quantum mechanics," 
in Doebner, H.-D., Scherer, W., and Schultz, C., (Eds.), GR0UP21. 
Physical Applications and Mathematical Aspects of Geometry, Groups, 
and Algebras, Goslar, Germany, June-July 1996 (World Scientific, Sin- 
gapore, 1997), Vol. 2, p. 764; "Transition probabilities between energy 
levels in the framework of the classical approach," invited lecture pre- 
sented at the Inaguration Conference of APCTP, Seoul, June 1996 (to 
appear in the Proceedings of the Conference, World Scientific, Singa- 
pore, 1998); "Classical description of quantum states and tomography," 



33 



in Han, D., Janszky, J., Kim, Y. S., and Man'ko, V. I., (Eds.), Fifth In- 
ternational Conference on Squeezed States and Uncertainty Relations, 
Balatonfiired, Hungary, May 1997 (NASA Conference Publication, God- 
dard Space Flight Center, Greenbelt, Maryland, 1998), Vol. NASA/CP- 
1998-206855, p. 523. 

[13] Man'ko, V. I., J. Russ. Laser Research (Plenum Press) 17, 579 (1996). 

[14] Man'ko, V. I., "Quantum mechanics and classical probability theory," in 
Gruber, B., and Ramek, M., (Eds.), Symmetries in Science IX (Plenum 
Press, New York, 1997), p. 215. 

[15] Man'ko, Olga, "Tomography of spin states and classical formulation of 
quantum mechanics," in Gruber, B., and Ramek, M., (Eds.), Symme- 
tries in Science X (Plenum Press, New York, 1998, to appear). 

[16] Mancini, S., Man'ko, V. I., and Tombesi, P., Quantum Semiclass. Opt. 
7, 615 (1995). 

[17] D'Ariano, G. M., Mancini, S., Man'ko, V. I., and Tombesi, P., Quantum 
Semiclass. Opt. 8, 1017 (1996). 

[18] Bertrand, J., and Bertrand, P., Found. Phys. 17, 397 (1987). 

[19] Vogel, K., and Risken, H., Phys. Rev. A 40, 2847 (1989). 

[20] Man'ko, V. I., and Safonov, S. S., Teor. Mat. Fiz. 112, 1172 (1997). 

[21] Man'ko, V. I., and Safonov, S. S., Phys. Atom. NucL 61, 585 (1998). 

[22] Dodonov, V. V., and Man'ko, V. I., Phys. Lett. A 229, 335 (1997). 

[23] Man'ko, V. I., and Man'ko, O. V., JETP 85, 430 (1997). 

[24] Leonhardt, U., Phys. Rev. A 53, 2998 (1996). 

[25] Dodonov, V. V., Malkin, I. A., and Man'ko, V. I., Physica 72, 597 

(1974). 



34 



[26] Man'ko, O. V., "Symplectic tomography of nonclassical states of a 
trapped ion," Preprint IC/96/39 (ICTP, Trieste, 1996); Los Alamos 
Report No. |quant-ph/9604018| ; J. Russ. Laser Research (Plenum Press) 
17, 439 (1996); "Symplectic tomography of Schrodinger cat states of a 
trapped ion," in Ferrero, M., and van der Merwe, A., (Eds.), Proceedings 
of the Second International Symposium on Fundamental Problems in 
Quantum Mechanics, Oviedo, Spain, July 1996 (Kluwer Academic Press, 
1997), p. 225; "Symplectic tomography of nonlinear Schrodinger cats of 
a trapped ion," in Han, D.,Janszky, J., Kim, Y. S., and Man'ko, V. I., 
(Eds.), Fifth International Conference on Squeezed States and Uncer- 
tainty Relations, Balatonfiired, Hungary, May 1997 (NASA Conference 
Publication, Goddard Space Flight Center, Greenbelt, Maryland, 1998), 
Vol. NASA/CP-1998-206855, p. 309. 

[27] de Matos Filho, R. L., and Vogel, W., Phys. Rev. Lett. 76, 608 (1996). 

[28] Nieto, M. M., Phys. Lett. A 219, 180 (1996). 

[29] Meekhof, D. M., Monroe, C, King, B. E., Itano, W. M., and Wineland, 
D., Phys. Rev. Lett. 76, 1796 (1996). 

[30] Haroche, S., Nuovo Cimento B 110, 545 (1995). 

[31] Man'ko, V. I., "Introduction to quantum optics," in: Castanos, O., 
Lopez-Pena, R., Hirsh, J. G., and Wolf, K.-B., (Eds.), Latin-American 
School of Physics, XXX ELAF, Group Theory and Its Applications, 
Mexico, July-August 1995, AIP Conference Proceedings 365 (AIP, New 
York, 1996), p. 337. 

[32] Man'ko, Olga, and Man'ko, V. I., J. Russ. Laser Research (Plenum 
Press) 18, 407 (1997). 



Fedele, R., and Man'ko, V. I., Phys. Rev. E 58, 992 (1998). 

[34] Cahill, K. E., and Glauber, R. J., Phys. Rev. 177, 1882 (1969). 

[35] Wiinsche, A., J. Mod. Opt. 44, 2293 (1997). 

[36] Man'ko, V. I. , and Mendes, R. V., "Noncommutative time-frequency 
tomography of analytic signals," LANL [physics/ 9712022] Data Analysis, 
Statistics, and Probability, IEEE Signal Processing (submitted, 1998). 

35 



[37] Smithey, D. T., Beck, M., Raymer, M. G., and Faridani, A., Phys. Rev. 
Lett. 70, 1244 (1993). 

[38] Dodonov, V. V., Man'ko, O. V., Man'ko, V. I., and Wiinsche, A., 
"Energy-sensitive and 'classical-like' distances between quantum states," 
Phys. Scr. (to appear, 1998). 

[39] Dodonov, V. V., and Man'ko, V. I., Invariants and Evolution of Non- 
stationary Quantum Systems, Proceedings of the Lebedev Physical In- 
stitute (Nova Science, New York, 1989), Vol. 183. 

[40] Malkin, I. A., and Man'ko, V. I., Dynamical Symmetries and Coherent 
States of Quantum Systems (Nauka, Moscow, 1979) [in Russian]. 

[41] Mancini, S., Man'ko, V. I., and Tombesi, P., Europhys. Lett. 37, 79 
(1997). 

[42] Mancini, S., Man'ko, V. I., and Tombesi, R, J. Mod. Opt. 44, 2281 
(1997). 

[43] Mancini, S., Man'ko, V. I., and Tombesi, P., Quantum Semiclass. Opt. 
9, 987 (1997). 

[44] Wheeler, J. A., and Zurek, W. H., (Eds.), Quantum Theory and Mea- 
surement (Princeton University Press, 1983). 

[45] Bell, J. S., Speakable and Unspeaskable in Quantum Mechanics (Cam- 
bridge University Press, 1987). 

[46] Wallentowitz, S., and Vogel, W., Phys. Rev. A 53, 4528 (1996). 

[47] Banaszek, K., and Wodkiewicz, K., Phys. Rev. Lett. 76, 4344 (1996). 

[48] Moyal, J. E., Proc. Cambridge Philos. Soc. 45, 99 (1949). 

[49] Landau, L. D., Z. Physik, 45, 430 (1927). 

[50] von Neumann, J., Mathematische Grundlagen der Quantenmechanik 
(Springer, Berlin, 1932). 

[51] Man'ko, O. V., Phys. Lett. A, 228, 29 (1997); SPIE (to appear, 1998). 

36 



[52] Man'ko, V. I., and Wiinsche, A., Quantum Semiclass. Opt. 9, 381 (1997). 

[53] Schrade, G., Man'ko, V. I., Schleich, W. P., and Glauber, R. J., Quan- 
tum Semiclass. Opt. 7, 307 (1995). 

[54] Schrodinger, E., Sitzungsber. Preuss. Acad. Wiss. 24, 296 (1930). 

[55] Malkin, I. A., and Man'ko, V. I., Phys. Lett. A 32, 243 (1970). 

[56] Malkin, I. A., Man'ko, V. I., and Trifonov, D. A., Phys. Rev. D 2, 1371 
(1970). 

[57] Man'ko, V. I., Rosa, L., and Vitale. R, Phys. Rev A 58, 3291 (1998). 

[58] Man'ko, V. I., Marmo, G., Sudarshan, E. C. G., and Zaccaria, 
P., "f-Oscillators," in Atakishiyev, N. M., Seligman, T. H., and 
Wolf, K.-B., (Eds.), Proceedings of the Fourth Wigner Sysmposium, 
Guadalajara, Mexico, July 1995 (World Scientific, Singapore, 1996), 
p. 421. 

[59] Man'ko, V. I., Marmo, G., Sudarshan, E. C. G., and Zaccaria, P., Phys. 
Scr. 55, 528 (1997). 

[60] de Matos Pilho, R. L., and Vogel, W., Phys. Rev. A 54, 4560 (1996). 

[61] Mancini, S., Man'ko, V. I., and Tombesi, P., J. Mod. Opt. 44, 2281 
(1997). 

[62] Mancini, S., Man'ko, V. I., and Tombesi, P., Quantum Semiclass. Opt. 
9, 987 (1997). 

[63] Andreev, V. A., and Man'ko, V. I., JETP 87,239 ( 1998). 

[64] Andreev, V. A., Man'ko, O. V., Man'ko, V. I., and Safonov, S. S., J. 
Russ. Laser Research (Plenum Press) 19, No. 4 (1998). 

Man'ko, V. I., Rosa, L., and Vitale, P., Phys. Lett. B 349, 328 (1998). 



37 



